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1.  MUNITION  EXPENDrPORE 

•—  — '  . .  . .  ' 

In  the  previous  progress  report  we  Indicated  that  the 

r 

Calder  Salzer  technique  for  inverting  a  Laplace  Transform 
iiiq>oses  IUnknown  conditions  on  the  moments  of  the  multi-munition 
dosage  d:\strlbutlon. 

or  course 4  as  Calder  i>olnts  out  a  direct  use  of  the 
Salzer  scleme  vrould  liiq;>ly  thtat  all  moments  are  infinite.—— 

Indeed  it  \/otad  imply  that  the  total  area  under  the  "distribution” 
is  infinite.  That  is  if  we  assume  that  the  Laplace  Transform, 

(1)  ,  ^(t)  =  /  e"^^f(D)dD  =»  2  ^  , 

0  •  lal  t^ 


then  the  djsage  distribution  is 


(2) 


2 

r»l^  (r-l)! 


Clearly  if  r(D)  is  given  by  (2)  then 


/  f(D)dD  =»  « 

0 


and  f(D)  cannot  bei  a  probability  density.  Actually  Calder 
uses  the  Laplace  Transform  of  the  cumulative  distribution  functloni, 
so  that  one  must  replace  f(D)  in  (l)  and  (2)  by  P(D).  As  he 
points  out 


P(D) 


k  b 

S  J£. - 

r«l  (r-l ) ! 


-  2  - 


(3) 


Implies  that  P(“)  =»  “  and  not  1,  as  vre  must  have  for  any 
cmnulatlve  distribution  functloa.  Differentiating  (3)  we  have 


f(D) 


k  b  0^-2 

Z  - - 

r=2  (r-2)! 


iTrom  v/hlch  X7e  drat/  the  same  conclusion.  Another  difficulty 
with  (2)  is  that  the  polynomial  goeq,  through  (0,0)^  Thus 
it  liaplles  zero  probability  of  obtaining  dosages  in  the  neigh'^ 
borhood  of  zero. 


Calder  handles  both  objections  as  follows.  To  get 

•  m 

arotind  the  non-zero  probability  of  zero  coverage  he  truncates 
the  dosage  distribution  at  the  lower  end.  To  get  around  the 
unboundedness  of  the  polynomial  at  D  »  «  he  Introduces  a 
scallxig  factor  (so  that  the  modified  P(D)  Is  bounded  by  1). 
Thusj  in  effect  he  uses  a  different  polynomial  fit  for  each 
D'^value.  It  Is  difficult  to  estimate  the  effect  that  this  mod^^ 
Iflcatlon  has  on  the  original  physical  problem.  Hov/ever,  It  is 
reasonable  to  asstaoe  that  a  family  of  polynomials  can  be  made 
to  fit  a  function  better  than  a  single  polynomial. 


In  addition^  although  the  scaling  technique  enables  one 
to  fit  an  P(D)  v/lth  the  property  that  P(D)  -*■  1  as  D  w  ^ 
one  would  probably  need  a  stronger  scaling  factor  to  guarantee 
that  all  moments  exist.  One  still  has 


3 


A  . 

llm  /  D^f(D)dD  »  00  , 

A  — »■  M  0 

so  that  this  method  really  lRq;>ll6s  that  the  moments  of  the 

* 

dosage  distribution  do  not  exist. 

It  should  also  be  pointed  out  that  the* theoretical 
Justification  for  fitting  f(D)  by  a  polynomial  5^  related  to 
the  Wfclerstrap  Approximation  Theorem  ifhich  states  that  any  con“ 
tlnuous  function  can  be  closely  approximated  by  a  polynomial  In 
D  for  real  D  satisfying  0  <  a  <  D  <  b  .  Notice  that  this 
theorem  holds  for  functions  over  finite  intervals ^  but  that  we 
are  dealing  with  Laplace  transfoxmis  which  are  Integrals  defined 
over  an  infinite  range.  We  can  then  show  that  by  using  the  un^ 
scaled  Laplace  Transform  method  It  is  not  possible  to  “fit”  the 
bounded  functions  by  a  polynomial.  To  rlgorlze  the  scaling  tech¬ 
nique  we  need^  a  corrollary  to  the  Welerstrap  Theorem^  v/hlch 
shov7S  that  a  sequence  of  polynomials  will  do  the  trick.  The 
Colder  technique  esqploys  a  sequence  of  polynomials ,  each  of 
which  fits  at  one  of  the  sequences  of  D-values. 


ii  In  the  last  progress  report  (po2)  we  wrote  the  Laplace  Transform 
in  terms  of  moments  of  the  distribution.  In  analogy  with  the 
moment  generating  function  \te  should  have  used 

<<(t)=  z  (-1)1—3^  . 

1=»0  1» 

The  (-1)^  v/as  Inadvertently  omitted. 
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Figure  1.  Polynomial  approximation  to  p(D) 

\}e  are  nov;  In  the  process  of  trying  a  less  sophisticated. 

?5 

technique  .than  scaling  to  allovf  for  the  fact  that  a  polynomial 
is  infinite  at  D  =*  «  .  It  consists  In  fitting  th&  polynomial 
betvreen  D  «=i  0  and  D  =  Dq  and  then  arbitrarily  approximating 
P(D)  by 

• 

P(D)  -  1  ,  D  >  Dq  . 

If  Dq  Is  large  enough  so  that  we  are  really  not  Interested  in 
the  probability  of  getting  dosages  above  Pq  and  if  our  approx°  * 
Imatlon  gives  P(D)  <1  D  <  ,  this  procedtire  may  be  satisfactory^ 
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2,  ATMOSPHERIC  DIProSIOM 

In  an  attempt  to  ol^aplify  the  dlfi^islon  model  treated 

In  those  reports  a  version  which  depends  on  dimensionless  para*- 

# 

meters  is  proposed^  As  a  by-product  of  this  ne\i  procedurt-^  we 

•  • 

are  enabled  to  obtain  a  family  of  depositions  by  raultlplying  by 
simple  expressions  the  result  of  integrations^  similar  to  (10) 
and  (11)  of  the  previous  progress  reports 


Our  deposition  of  agent,  K  as  a  function  of  distance 

r 

doimwlnd,  x  v/as  given  by 

^2 

(5)  K(x)  =  A  /  h(n)t(x,0  dr/ 

«1 

.v;here 


(6)  h{^25^i,,t)  =  —  exp[- 

^2?r)2  t  at® 


t(x,n) 


— (I)'""' 

r(i+p;  * 


“f/x 


(8)  p  =>  cO 


®  \Ie  nov;  make  tv/o  changes.  In  place  of  x  we  consider 
the  dimensionless  parameter  x/f  ,  For  any  value  of  f  it  is 
then  a  simple  matter  to  treat  K  as  a  function  of  x.  Secondly 
the  integration  Is  performed  with  p,  rather  than  Cl  as  the 
variable  of  Integration.  Thus 


~  6  “ 


(9)  ~  /  h(pj(p,t)  v(r,p)dp  t=»  M(r59,t)  ,  say 

Ac  0 

vihere 

(10)  r  =  x/f 


(11) 

(12) 

(13) 

« 

The  meteorological  parameters  are  toeorporated  in  A,  c 

and  f  =  Q/(Ac).  Thus,  each  curve  M  can  be  easily  transformed 

by  (10),  (11)  and  (12)  to  a  set  of  curves,  K,  for  specific 

* 

meteorological  conditions  and  particle  size  distributions.  The 

conditions  treated  in  previous  reports  are  given  by  A  =  .1674  ~ 

m*  see 

and  the  values  in  Table  1.  It  should  be  noted  that  in  all 
these  cases  9  ,  the  mean  of  the  “normal”  distribution  in  (11) 
v/as  negative.  Since  9  «  rather  thai  \i  seems  to  be  the  per¬ 
tinent  “mean”  to  consider  in  evaluating  K  it  was  decided  to 
compute  some  values  for  9=0  and  9=1  ,  where  we  expect  diffusion 

«  7 


l^(P5®^t)  =  ^  expE 

t 

-  3- 

"ySF  t 

.  9  ss  p,  «  jgnc 

v(r,p)  =  [r(Hi?)  r^'^P  e^  3 


to  play  a  larger  role.  For  thlo  reason  M  ,  of  (9)  Is  being 
evaluated  on  the  IBM  650  for  the  parameters  of  Table  2  by 
Mr.  Daniel  H.  Sallcoff,  who  performed  the  previous  calculations. 


f 

• 

c“^ 

j6x)C 

t  . 

9 

51^ 

4.06 

.08603 

1.7595 

.53 

f  2.30 

514 

3.26  ‘ 

.08603 

1.7595 

.57 

•^-1.50 

1028 

4.06 

,17206 

2.453 

.53 

•h  1,61 

1028 

3.26 

,17206 

2.453 

.57 

-h  ,81 

Table  1,  Parameters  In  previous  deposition  calculations . 


9 

t 

Pi 

P2 

0 

.3 

,40 

2.5 

0 

.55 

,19 

5.2 

0 

1.0 

.05 

20,0 

1 

.55 

CQO- 

jrC't.t- 

.071 

l4-7i3^ 

/.  9 

Table  2,  Values  of  parameters  being  used  In.  computing  M  as 
given  in  Equation  (l4),  where  **  ®  » 


O^ls  nev;  prooedure  outs  out  the  nsoesslty  for  evaluating 
L  ,  the  Integral  iflth  p  and  f  (or  r)  replaced  by  2pg  2f 
(or  If  one  actually  evaluated  (9)  exactly «  the  change  In 

these  parameters  would  merely  change  the  Integrating  grid  aod 
give  different  values  of  distance  dotmtilnd.  However «  we,  of 
course,  replace  (9)  by 

(14)  /  b(P*<P#t)  v(r,p)dp 

Pi 

where  and  p^  are  chosen  so  that  the  contribution  to  M 

of  values  \iheve  Pg  and  below  p  ,  are  negligible.  Replacing 
p  by  2p  v/ould  also  change  the  limits  of  Integration,  thereby 
omitting  liiqportant  values  of  p  and  replacing  them  by  values 
which  make  virtually  no  contribution  to  M  , 

An  outline  of  a  nev;  stochastic  model  of ,  turbulent 

diffusion  Is  given  In  an  appendix  to  this  report.  Of  course,  ' 

the  details  eire  yet  to  be  prograosed  for  coiig>utatlon  on  a  digital 

counter,  and  these  results  would  have  to  be  compared  vflth  field 

test  data  to  see  how  well  the  model  fits. 

*■ 

a 

V 

Leon  H.  Herbach 
Project  Director 
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my  195:1. 

A  stochastic  model  of  tiirbulent  diffusion 

The  purpose  of  this  memorandum  is  to  formulate  a  nevr  model 

of  turbulent  diffusion  by  a  simulation  of  the  physical  process,, 

The  model  v;ill  describe  the  process  dtu’-lng'a  fixed  Instant  of 
•  It  may  be  described 

time.  In  mox^  picturesque  language, as  being  a  "snapshot"  of 

the  physical  process c  It  is  hoped  in  later  reports  to  describe 

a  continuous  model  for  the  same  physical  phenomena. 

It  is  assiuaed  that  an  instantaneous  point  source  located 

» 

at  X  Xq  •-  0  ,  y  =>  =»  0  originates  as  a  "mass"  Qq  x^hlch 

my  be  described  mcithematlcally  by  the  delta  function,  or  unit 
impulse  .function,  Qq  3(x,y)  v/he-re  6(x,y)  is  defDJied  by 

<5{x,y)  =0  if  X  0  or  y  ^  0 
(1)  '^(x>y)  -  if  X  =  0  and  y  =  0 

+03 

Qq  ''  ^  *(2c,y)  =  Qq  ’ 

Our  .next  ascTimption  is  that  the  mass  x/111  diffuse  according 
to  a  "rar.dora  xjalk"  process  vihose  mechanism  is  described  below,, 

1/e  shall  be  pax-'ticularly  interested  in  studying  the  correlation 
functions  of  the  stochastic  process  involved.  Tiose  correlation 
functions  will  be  used  to  describe  the  rate  of  diffusion,  and* 
the-  general  features  of  the  process  mechanism.  The  con’c-lation 
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function  as  generated  vjlll  consider  the  relations  hetv/een  fixed 
discrete  points . 

It  is  assumed  that  the  "mass"  undergoes  a  splitting  vjhich 
gives  rise  to  tuo  "daughter”  masses,  liach  of  these  daughters 
iias  a  total  ’'mass”  Qq/2  and  is  cente7.’ed  at  the  points 

12)  =*25;^  +  6^^^-  5  =  ^0  1=1,2 

v;here  and  are  normally  and  independently  distributed 

v:ith  meauis  zero  and  variances  a£  . 

‘fhe  "splitting”  process  continues  in  this  fashion  except; 
that  the  variances  of  the  nev/  e’s  and  6«b  are  o|  , 

1  =  2, .3, 4, 5,  vjith  cr^  <  <  a|  <  a|  <  a|  ,  until  v^e  arrive 

at  the  points  v/here  =  1,2,..., 32.)  The 

process  is  illustrated  in  Figure  1.  That  is,  the  process  con¬ 
tinues  until  v;e  have  generated  2^  or  32  points. 

The  effect  of  v/ind  can  be  talcen  into  account  by  this  model. 

If  a  constant  v/lnd  is  blov/lng  from  some  direction,  its  effect 
on  the  distribution  of  the  sample  points  can  be  alloi^ed  for 
by  ass’imlng  tiuit  the  mean  of  the  first  points  ,  6^  ( 1=1,2) 

are  non-zero.  For  example  if  the  v/lnd  is  from  the  SW  we  may 
talcc*  E(€^)  =  1  ,  2(^i)  =*  i  » 

At  each  of  the  32  x>olnts  which  have  been  generated  we 
now  conceive  of  the  mass,  as  being  distributed  according 
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\ 


/■ 


(4‘^>4‘^h 


I  (4^K  yi^h-^/  '(x^>.y^>),i  _ 

!  '  f 

‘  t 

X  »/>),/- - ■'  -  ..  /'■ 


(x(^),y(l)),; 


(xi“>,4^‘‘>) 


(4®^yF>) 

J 


Figure  1.  The  Splitting  Process 
Oeneratlng  the  Random. 
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to  a  squared  exponential.  This  Is  tbe  llioltlng  result  If  we 

•  * 

continue  the  splitting  ad  infinitum  with  small  variances, 
^usj  the  density  at  any  point  (u«v)  in  spaoe  is 


+  (v 


where  the  points  are  generated  In  the  following  manner; 


i^Q*yo) 

-  (x("> 

.  yp>) 

.  (xp^ 

.  ypb 

(x<^)  , 

— 

(4^). 

y(^h  , 

(*2®^  . 

yi^h 

ypb  - 

(x<">  . 

y</>)  , 

(4"). 

yi^^b 

{x(^)  . 

4^b-^ 

(x|^^  , 

y^^>)  . 

(4®^  - 

y^^b 

** 

(4*’  - 

yi*>)  . 

(xf’  . 

y^b 

(x|^>  . 

- 

(45)  , 

yi^h  . 

{4®>. 

4^b 

(4“^ 

yi''>)  - 

(4">- 

y^’)  . 

(4®>- 

y'fb 

The  reason  for  using  a  squared  exponential  as  the  limiting 
distribution  for  the  "splitting"  of  steps  6  to  Infinity  - 

the  ..’fijul';  ing  dlstrlbtxtlon  is  aingulaTf  vz  folX  /'ra; 

aiiount  iti  any  region  after  an  Infinite  mutf.ci’  of 
3.6  given  bj-  the  circular  rorjii'il  integral.  Us  are  3.>:j 
tl<£  Vii/ce  after  step  five;,  by  their  a.verrigec.  It 
io  clt.:-av^  from,  tbs  v^-ay^  tL®  variances  of  subsequent  stepr  de'vreasf 
that  the  vai*Iab.ility  in  these  .latter  steps  is  neglig.lble  cmiiarf'd 

with  the-  «^:.xf*lier  stops.  Actual  calculations  will  Indicate  whetbei 

-  13  - 
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a  good  approxliaatlon  will  be  obtalziad  If  this  ''averaging*'  is  done 
after  step  5#  6,  or  some  other  nijmber. 

It  should  also  be  noted  that  we  have  still  not  taken 

time  into  account «  since  all  stages  In  the  splitting  in^ocess 

* 

are  assumed  to  take  place  instantaneously.  Time  is  considered 
by  replacing  in  (3)  by  o*f(t)  ,  vihere  f(t)  is  an  in« 
creasing  function  of  time.  The  first  function  to  be  considered 
will  be  f(t)  -  t. 

e 

Noiv  we  omy  consider  the  correlation  structure.  It 
may  be  noted  that  the  o^s  idiioh  nust  be  used  here  will  have 
to  be  obtained  experimentally.  The  correlation  structure  for 
the  first  two  steps  is  given  as 


(9) 


og  0 

^  0  og 


’oh 


In  considering  the  diagram  shown  in  Figure  2  it  is  noted  that 

if  the  correlation  structure  is  given  in  the  following  mauuier 

of  describing  the  splitting  process  we  have 
•• 

(10)  X-Xq-0  ,  y-y^-o  With  Oq 

as  the  splitting  continues  we  obtain  at  the  first  step  the 
following  o*s  ,  whexHi  Xq  »  flz^t  mass  and  x^  ,  the  second* 
mass,  considering  all  possible  combinations 
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IS/f 


‘■00 


:n) 


^01 


‘•10 

0 

0 


‘‘11 

0 

0 


and  at  the  next  step  we  obtain  the  following: 


*000 

*001 

*010 

*011 

*100 

+  o|  +  o| 

'o*®! 

»s 

0 

(12) 

og  +  of  +  »| 

»s 

»s 

0 

og  +  o«  +  o|  • 

og  +  a® 

0 

ff®  +  CT®  +  O® 
0  1  2 


Now  consider  the  correlations  of  these  realizations: 


(13)  At 


xiihere 


Rq  -  I 


'i  ^ 


’Bh 


*  t 


(Identity  matrix) 


1 

3 

1 


16 


or 


13.8 


*  ^.0  t  ^  ^ 


(14) 


of  +  o»  Eo  •  I9 

^  [  0  _2  2  1  2  -2, 


+  Og  Ig  •  Ig 


It  Will  be  noted  that  the  covariance  natrlx  at  any  given 

*  •  it  » 

realization  of  the  prooeas  gives  the  history  of  the  mass"  up 

to  this  point «  and  essentially  describes  the  fact  that  a  hlg^ 

correlation  will  be  found  for  points  which  lie  In  a  close  re> 

latlon  physically  to  each  other.  Vihlle  they  have  a  low  oorz^» 

» 

latlon  with  the  "aass"  points  which  lie  at  some  distance  from 
each  other.  This  Is  noted  In  Figure  2  and  seems  to  agree  with 
the  physical  theory. 

To  end  the  mattdr  of  notation  we  will  designate  the  points 
at  any  given  realization  In  their  natural  order.  At  the  fifth 
step  we  will  nisBber  the  z*eallzatlons  as  ...^x^g  . 

Of  course  •  the  details  of  this  model  still  have  to  be  worlsed 
out.  Among  the  problems  to  be  Investigated  are:  (l)  How  can 
one  obtain  the  mean  variance  and  the  distribution  of  dosage  at 
a  point  or  set  of  points?  (11)  Uhat  is  the  limiting  distribution 
of  dosage  at  a  point?  If  we  drop  squared  exponentials  uniformly 
over  a  large  target  the  resulting  distribution  is  c.d.  We  ti:ill 
probably  be  involved  with  a  distribution  which  Is  the  sum  of 
c.d®3. 
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